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Abstract 

Consider the one-parameter generalizations of the logarithmic and exponential functions which are obtained from the integration 
of non-symmetrical hyperboles. These generalizations coincide to the one obtained in the context of non-extensive thermostatistics. 
We show that these functions are suitable to describe and unify the great majority of continuous growth models, which we briefly 
^ review. Physical interpretation to the generalization function parameter is given for the Richards' model, which has an underlying 
microscopic model to justify it. 
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1. Introduction 



X 



The convenience of generalizing the logarithmic and 
exponential functions has attracted the attention of re- 
searchers in the last years. One-parameter generalizations 
of logarithmic and exponential functions have been pro- 
posed in the several contexts |1I2I3I4I5I6I7| . More sophis- 
ticated generalizations have also been proposed such as a 
two and three-parameter |8l9j . which includes all the oth- 
ers generalized logarithms as special cases. Examples of 
the convenience of these generalizations have been seen in 
different fields |10I11I12I13I14I15I16I17I1"8] . 

From an initial population size no, at an initial time to, 
growth models are used to predict the size n(t) of popu- 
lation at an arbitrary time t > to- The most important 
parameters are the intrinsic growth rate r and carrying ca- 
pacity = n(t — > 00). Nevertheless frequently additional 
parameters are needed to adjust the model to a given ex- 
perimental data. Up to date, a comprehensive description 



of continuous growth models is given by the Tsoularis and 
Wallace model [TH] , which generalizes the logistic- like mod- 
els [20l2ll22l23l24l25l2B] . 

The main objective of this paper is to show that the 
one-parameter generalization of the logarithmic and expo- 
nential functions proposed [T] is suitable to describe con- 
tinuous growth models. This paper is structured as fol- 
lows. In Sec. [5J we present the one-parameter generaliza- 
tion of the logarithmic and exponential functions and a 
brief review of growth models. In Sec. [3] we consider a 
simple type of Bernoulli's equation, which solution is the 
Zipf-Mandelbrot function. From the gradual generalization 
of the differential equations, their solutions retrieves the 
Richards' growth model [55] , Marusic and Bajzer model 26J 
and finally Tsoularis and Wallace model [19]. Since the 
Richards' model has an underlying microscopic model |27j . 
one is able to interpret physically the logaritmic function 
generalization parameter. Final remarks are presented in 
Sec. II 
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2. Brief review 

To facilitate the connection between the one-parameter 
logarithmic and exponential generalizations and growth 
models, we briefly review these subjects in the following. We 
first present the simple geometrical arguments [28] which 
we have used to obtain a one-parameter generalization of 
the logarithmic and exponential functions. This natural ge- 
ometrical generalization coincides exactly to the one pro- 
posed by Tsallis [1] in the context of non-extensive ther- 
modynamics [29] . Instead of calling the extra parameter as 
q, like in the precursor papers, we have used q = 1 — q, 
which is more convenient to stress the symmetry properties 
of these functions. Next, we present some general continu- 
ous growth models and show how to retrieve the standard 
ones with the appropriate set of parameters. 

2.1. The q-logarithmic, and the q- exponential functions 

The ^-logarithm function ln^(a;) is defined as the value 
of the area underneath the non-symmetrical hyperbole 
fg(t) = l/i 1 ^ in the interval t £ [l,x], with x > 0: 
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q'^q 



x q - 1 
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(1) 



For any q, the area is negative for < x < 1, it vanishes for 
x = 1 [hiq(l) = 0], and it is positive for x > 1. This func- 
tion is not the logarithmic function in the basis q [namely 
log^x)] but a one-parameter generalization of the natural 
logarithmic function definition. 

The ^-exponential function eq(x) is defined as the t- value, 
in such a way that the area underneath fq(t) — 1/t 1- *, in 
the interval t £ [l,e^(x)], is x. This is the inverse of the 
g-logarithm function eq[hig(x)} = x = \n.q[eq(x)\ and it is 
given by: 



eq(x) = Jimjl + q'x]+^ 



(2) 



Here the the operator [a]+ = max(a,0) is necessary since 
eq(x) is not real if qx < —1. This is a non-negative function 
Sq{x) > and x = is a special point since ej(0) = 1, 
independently of the q value. 

2.2. Continuous growth models 

Population dynamics deals with the evolution of species 
interacting with themselves and the environment. The 
one-species (logistic-like) models are well described in 
terms of growth models. Essentially, growth models con- 
nect the relative growth rate dn(t)/[n(t)dt] = dhin(t)/dt 
of a population, with n(t) individuals of a given species 
at time t, to an induced saturation function G r/llao (n): 
dlnn/dt = G r . noo (n), where — n(oo) is the carry- 
ing capacity (environment saturation level), which corre- 
sponds to the asymptotic value of the stationary solution 



[dlnn(t)/dt = G r . noo (n) = 0] and r is the intrinsic growth 
rate (individual contribution). To compare populations in 
different environments, a better variable to work with is 
the ratio between the population size and the carrying ca- 
pacity (final population size): p(t) = n(£)/noo. In this new 
variable, one has: 



d hip 
dt 



G r (p) , 



(3) 



where we notice that now, the induced saturation function 
depends only on r, which we write G r cx r. Thus, the envi- 
ronmental contribution is taken into account into the vari- 
able itself. 

In the following we consider one-species (logistic-like) 
models. The Malthus' (exponential) model is the simplest 
one [G r {p) — r L which is extended basically to two other 
models: 

- the Verhulst (logistic) model [G r (p) = r(l — p)], whose 
discretization leads to the logistic map (prototype exam- 
ple in dynamical system) and 

- the Gompertz's model [G K (p) = — nhip], which is also 
used to calculate life insurances taxes [20|21j . 

It is interesting to observe that in Ref. 30, p. 83], it 
is shown that the Verhuslt and Gompertz models can be 
generalized by a simple function proposed by Montroll 
and Badger [31] . which in fact is exactly the Richards' 
model 1221 



G K (p) 



l-p*(t) 



(4) 



where k = rq. When q = 1 , one obtains the Verhust model 
and when q — > 0, one retrieves the Gompertz's model. It 
is interesting to point out that k and r are not exactly the 
same intrinsic growth parameter, although they are linearly 
related. 

A more general description of growth models has been 
given recently by the Tsoularis and Wallace model pjU 
p. 22]. They considered the following induced saturation 
function: 



G R [p) = Kp a - l {t) 



i- P Ht) 



(5) 



The formal solution of Eq. ([3]) with 
the above induced saturation function is given by: 



where k = ripn 1 ^ 
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, 1 - 7 = «t , (6) 



where B x (a, b) — dt i a_1 (l — t) b ^ 1 is the incomplete beta 
function and two sets of values are acceptable: 

(i) a < 1, q > and 7 < 1, 

(ii) a > 1, q < and 7 < 1. 

Consider the special cases of the above model. 
- For 7 = 0, one has Zipf-Mandelbrot-like growth models, 
which has two special cases: 

• exponential, for a = 1 and 

• power-law, for a > 1. 



For 7 = 1, one retrieves the Marusic and Bajzer 
model 26J, which generalizes other growth models: 

• if q = 1 — a, the generalized von Bertanffy's model or 

• if a = 1 and q > — 1, the Richards' model. 
For other 7 values, one obtains: 

• For a — 1 and <7 — > 0, the hyper-Gompertz's model 
G K (p) = -Adn^(t). 

• For g = 1, one obtains the Blumberg's growth 
model [24] used to model organ size evolution. The 
Blumberg's model reduces to the hyperbolic model for 
regenerative growth, if 7 = 1 + 1/iV and a = 1 — 1/N 
(seeRef. [23]). 

■ If a = 1+9(7— 1) , Tsoularis and Wallace model reduces 
to the Turner et al. model, which surprisingly admits 
analytical solution [25] . 



3.2. Richards-like growth models 

Eq. ([7]) can be generalized to take care of other interesting 
situations. The following differential equation: 



dj/Qr) 
dx 



my(x) + (k — rnjy 1 q (x) , 



(9) 



has been used in growth models [26j and other con- 
texts |32I33| . Factoring my{x) and calling y = [m/(m — 
k)] 1/q y(x), Eq. © is written as: 



dlny (x) 
dx 



= mq In, y 1 (x) , 



(10) 



which leads to the Richard's model as we will show. 



3. (/-generalized functions and growth models 

In what follows we will first deal with a particular type 
of Bernoulli's equation [d/dx + p(x)]y(x) — q(x)y 1 ~ q (x), 
with p(x) = and q{x) = 1. We show that this equa- 
tion (kinetic equation of arbitrary order) is related to the 
hyper-Gompertz model. When a new linear term is added 
to this equation, it represents the Richard's growth models 
and in this case we have been able to give a physical in- 
terpretation to q. If this new added term is non-linear, we 
obtain an equation, which represents the Marusic-Bajzer 
models. The effort term can be included considering the q- 
logarithm function raised to a power, which gives rises to 
the Tsoularis and Wallace model. 



3.1. Zipf- Mandelbrot model 

The derivative of the ^-exponential function with re- 
spect to x is: deq(kx)/dx — k [e^fex)] 1 "*, so that the q- 
exponcntial function is the solution of the following non- 
linear first order differential equation [32133135] : 



dj/Qc) 
dx 



= k y x -~«{x) 



(7) 



where k is a constant. In Chemistry it is known as kinetic 
equation of order l — q. This equation is not directly related 
to any growth model, but if one thinks that y = hip, then 
one has the hyper-Gompertz model. 

The solution of Eq. ([7]) can be found rewriting it as [2JJ : 
qy q ~ l dy/dx = qk =>■ dy''(x)/dx = qk, which is linear in the 
variable v(x) = y q (x) and can be solved: v(x) — qkx + vq, 
where the "initial" condition is vq — v(0) = y$ with yo = 
y(0). Thus, one writes: 



with Xq — y\jk. 



1/9 



(8) 



3.2.1. Richards 'growth model 

In Eq. (|10p . calling t the independent variable x, p = l/y 
and mq = — k, one rewrites it as: 



d\np(t) 



dt 

which solution is 
Pit) ' 



K lUqp(t) 



(11) 



-l\„-reti 



e^q[\n-q{p )e- Kt } 



e_^[-ln 5 (> )e 



(12) 



This is the Richards' growth model [22] . 

The discretization of Eq. [TTJ leads to a generalization of 
the logistic map different from the previously considered 
one (Feigenbaum's map). An initial study of this recursive 
equation has been presented in Ref . [33] . 

3.2.1.1. Particular Cases. Some particular cases can of 
Eqs. (11]) and JTJ]) be retrieved: 

- If q = 1, one recovers the Verhuslt's model (logistic equa- 
tion): 



dlnp 
dt 



= -k hii(p) = k (1 -p) , 



which has the following solution: 



p(t)=e- 1 [(l-po 1 )e- 



1 + (Po 1 - l) f 



(13) 



(14) 



where pa = rio/ n oc with no = n(0) being the initial 
population size and n(oo) = is the carrying capacity. 

• If k > 0, then for t ^> 1/k, p(oo) = 1, i.e. the carrying 
capacity is an attractor. 

• If k < 0, the carrying capacity tIqq becomes a repulsor 
and the solutions p(t) = l/[l+bo 1 -l)e |re|t ] = e_i[(l- 
Po 1 )e' K ''] ) either converge, for n < (i. e. p < 1) 
or diverge, otherwise. This divergence means that the 
population becomes non-limited for times greater than 
t* = ln(l — Pq 1 )/ k, which is solution of 1 + (p^ 1 — 
l)e Kt =0 and depends on the initial value Pq. 



As rioo — > oo, then p — > 0, and one has: 
dlnp 



(15) 



which is the Malthus' model. The solution of Eq. (fT5")) 
is: 

p(t) = p a e Kt . (16) 
If q = 0, one retrieves the Gompertz's model: 
dlnp 



dt 

which solution is: 
p(*)=Po~" ■ 



(17) 



(18) 



If q = — 1, one has the monomolecular or Mitscherlich 
growth model [HI p. 33]: 



dp(t) 



= fc[l-p(t)] , 



dt 
so that 

p(t) = l-(l-po)e- 



-it 



(19) 



(20) 



Observe that Mitscherlich model [Eq. (|19pj. which 
uses dp/dt differs fundamentally from Verhulst model 
[Eq. (|13[)]. which uses d\np/dt. 

3.2.1.2. Microscopic Model. The competition between 
cell drive to replicate and inhibitory interactions, that are 
modeled by long range interaction among the cell, furnishes 
an interisting microscopic mechanism to obtain Richards' 
model |27j . The long range interaction is dependend on 
the distance r between two cells as a power law r 7 and 
the cells have a fractal structure characterized by a fractal 
dimension Df. Using Eq. [1] one can write their Eq. (7) of 
Ref. [27] as: dlnn(i)/di = (G) - Juh^[Dfn(t)/u], where 
w is a constant related to geometry of the problem and 
q = 1 — j/Df. Here the parameter q acquires a physical 
meaning related to the interaction range 7 and fractal 
dimension of the celular structure Df. This physical in- 
terpretation of q has only been possible due to Richards' 
model underlying microscopic description. 

Another point of view is to obtain particular cases from 
the microscopic model of Ref. [27] . If all the cells interact 
with the same intensity for a fixed Dt ^ 0, then 7 = 
and q = 1, so that Eqs. (fTTjl and (fT2"|) lead to the Verhulst 
model. Another way to obtain the Verhulst model is to fix 
the cell interaction range 7^0 and consider large fractal 
dimensionality Df — > 00. Both cases illustrates that the 
Verhulst model is a mean field model. If cell interaction 
range matches with the cell fractal structure then 7 = Df 
and q = 0, so that Eqs. (|11|) and (fT2"|) lead to the Gom- 
pertz's model. If 7 > Df, then q < and Eq. (fT2|) reads: 
p(t) = 1 + hi|,j|(po)e~ Kt , for nt ^> 1, which leads to expo- 
nential behavior. In particular if 7 = 2Df, one retrieves the 
Mitscherlich model [Eq. {20])]. 



3.2.2. Richard- Schaefer's model 

If individuals can be inserted to or removed from a given 
population at fixed time intervals at a rate e, which is called 
the effort rate, then one has the Schaefer's model [20|21j : 

d lnp 



dt 



-k hig(p) - e 



(21) 



The solution of Eq. (J2TJ> is: 

e 5 (e) 



Pe, K (t) = 



eg{ln 5 [e g -(e)po 1 ]e-(«-^)'} 
= e- q {e) e-~ q {- We^K^e"'"-**} , (22) 
so that when e = 0, one reobtains Eq. (fT2"|) . 

3.3. Marusic and Bajzer's growth model 

As proposed by M. Marusic and Z. Bajzer [55] the equa- 
tion 



—j^- = ay (x) + 6y p (x) 



(23) 



with a, b,a,(3 € 5i, furnishes the basis of the models of one- 
species dynamical population, which generalizes traditional 
growth models. 

We remark that Eq. (|23p with a > and 6 < was first 
studied in Ref. [36], which deals with quatitative laws in 
metabolism and growth. A more fundamental derivation of 
Eq. (|23|) has been furnished in Ref. [37 38 in the allometric 
morphogenesis (form or patern generation) context. Maru- 
sic and Bajzer have considered arbitrary real parameters 
and solved in Eq. ([23]) . 

Further, we notice that one can rewrite Eq. ([23]) as: 
dy(x) /dx = ay a {x) + by a (x) (J3 - a) [/" Q {x)-l + l]/{§- 
a) = ay a (x) + by a {x) lnp- a (y), with a = a + b and b = 
b(/3— a). In this way Eq. (2) in the Marusic and Bajzer pa- 
per |26] is generalized and retrived for the particular case 
/3-> a. 

Calling a = m, b = k — m, a = 1 + q' and /3 = 1 + q' — q 
in Eq. (|23]) . one obtains: 



dy(x) 
da; 



my 



1+ «'{x) + (k-m)y 1 



\x) 



(24) 



Taking q' = 0, one recovers Eq. ([9]). Factoring my 1+ ^ (x) 
and calling y — [m/(m — k)] 1 ^y(x), Eq. p4]) is written as: 



dx 



mq 



q'/q 



(25) 



Calling t the independent variable x, p = l/y and mq[(k 
m)/mf '* = k, one rewrites Eq. (|2"B"]) as: 



d ln^> p(t) 
dt 



-K ln.gp(t) 



(26) 



This equation has analytical solutions only for some par- 
ticular values of q' . For q' = 0, one recovers the Richards' 
model [Eq. For q 7^ 0, other common growth models 

arc obtained as limiting cases, as we see below. 



A 



3.3.1. Generalized von Bertalanffy's model 

For q' = q, one has the generalized von Bertalanffy 
model [19 , p. 47], and Eq. (j26|) has the following solution: 



p(t) = eq[liig(p )e 



(27) 



If q = 1/3, then q 1 = q = 1/3 and one recovers the special- 
ized von Bertalanffy model [36] . 

3.3.2. Smith's model 

The Smith's model [191 p. 35] replaces the term 1 — p in 
Verhulst model by 1 — (p + adp/dt). The resulting equation 



dlnp(t) 
dt 



l-p(t) 



1 + rap(t) 



(28) 



which cannot be obtained immediately from Eq. ([2"6"[) , but 
from an approximation of it by setting q = 1 and q' = 
1 - 0.473. 

3.4. Tsoularis and Wallace's growth model 

The growth model proposed by A. Tsoularis and J. Wal- 
lace [19] [given by Eq. (5])] can be obtained generalizing 
Eq. lHJJ} to: 



dhag p(t) 
eft 



K [- lllgKt)]" 



(29) 



which solution is given by Eq. ([6]), with a = 1 — q'. For 
7 = 1, one retrieves the Marusic and Bajzer growth model 
[Eq. (|26p ]. For 7/1, other common growth models are 
retrieved as limiting cases. 

The parameter 7 has the effect of an effort rate. This 
can be seen rewriting Eq. (|29[) in terms of a supplementary 
g-logarithmic function: 



dlligi p(t) 

dt 



= ft ln 7 [— hiqp(t)] — e 



(30) 



where k = and e — —k. For comparison, see Eq. ([2T|) . 

3.4.1. Hyper- Gompertz's model (or Exponential 
polynomial) 

For <f = q = in Eq. ([2"9"[) , one finds the hyper- 
Gompertz's model 



dlnp 
dt 



= k [— hip]' 



(31) 



Calling y = hip, this equation becomes Eq. ([7]), whose 
solution is given by Eq. ©. In the original variable, the 
solution of Eq. ([3~lj) is: 

p(t) =P0 1_o [_Kt/(lnP0)1_71 = e HPo)e!- -,[~Kt/ (In pa) 1 —* ] ^ ^2) 

where we point out the presence of a generalized exponen- 
tial function in the argument of an ordinary one. 



3.4.2. Blumberg's model 

For q = 1 in Eq. ([29|) , one retrieves the Blumberg 
model [24], which has been proposed to model organ size 
evolution: 



which can be rewritten as: 

P(t) 

<ip p 9 _1 (1 — p) ' = nt 



(33) 



(34) 



If q' > and 7 < 1, one has a closed analytical solution [as 
in Eq. ©, a = 1 - <f]: 



B p (t){q', 1 - 7) - -Bpo^ 1 - 7) = Kt 



(35) 



Notice that in Eq. (33]) for: 
-7 = 0, one has the power law induced saturation function 

G K (p) = Kp 1 ^^ and 
-7=1 + 1/N and q = 1 — 1/iV, one has the hyperbolic 

model for regenerative growth (see Ref. [2"3"]). 

3.4.3. Turner et al. model (Generic growth) 

For g' = (7(7 — 1) in Eq. ([2"9"| . one retrieves the Turner 
model 125 : 



which has the following analytical solution: 
P(t)=e-q {hi- q (po) ei_ 7 [-Kt(hi q Po 1 ) 7 " 1 ] } 



(36) 



(37) 



The Turner et al. model can be solved analytically only 
because of the judicious choice of q'. 

3.4.4. Retrieving growth models 

As one can see, all continuous growth models considered 
in this paper are generalized by Eq. ([2"9"]) . To recover those 
models, take appropriate values of q' ', q, 7 and k shown 
in Tab. [1] Some particular cases have analytical solution 
and they are shown in Fig. [1] It can be observed that for 
the Richard-Schaefer model, the asymptotic value is Poo = 
e q ( e ) = Vol = 0.89 . . ., due to the effort e = -0.1 rate, 
instead of poo = 1. 

4. Conclusion 

The main objective of this paper is to show that the one- 
parameter generalization of the logarithmic and exponen- 
tial function obtained from non-extensive statistical physics 
is suitable to describe a large class of continuous population 
dynamics (growth) model. On one hand, the g-logarithmic 
function appears in the differential equations that define 
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dlnp/dt = r(l — p) 
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dlnp/dt = —nlnp 


hyper-Gompertz 








* 
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din p/dt = /c(-lnp)7 


Richards 
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dlnp/dt = r(l -p«) 


Tsoularis and Wallace 






* 
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d Inqi p/dt = k(— Iriqp) 1 


Marusic and Bajzer 




* 
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d In^i p/dt = — k Inq p 


Mitscherlich (monomolecular) 


1 


1 
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lni dp/dt = dp/dt = k(1 — p) 


Blumberg 
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* 
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d \riqi p/dt = k(1 — p) 7 


Turner et al. 


1 + 9(1-7) 


* 


* 
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dlnp/dt = Kp^-^i-lnqpP 


Specialized von Bertalanfly 


1/3 


1/3 


1 


* 


d\n 1/3 p/dt = -Kln 1/3 p 


Generalized von Bertalanfly 


g 


* 


1 


* 


d\nqp/dt = —Klnqp 


Smith 


1 - 0.473 


1 


1 


* 


approximation 



Table 1 

Growth models which can be obtained from Eq. J29j. The asterisk means that the considered parameter can take an arbitrary value within 
its domain. 
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Fig. 1. Models which present analytical solution from Eq. <|29j l. The 
curves have been obtained with the following values: po = 0.001, 
a = 0.1 and e = —0.1. For parameters of arbitrary value in Tab. [l] 
we have taken: k = 1, q = 2, and 7 = 1.5. 

the models. On the other hand, the g-exponential func- 
tion appears in the solutions of those differential equations, 
when it is analytically available. In the presentation we have 
gradually increased the number of parameters of the differ- 
ential equation to establish connection with growth mod- 
els. Eq. (j2"5|) represents the most general model, the one of 
Tsoularis and Wallace. It has a simple form when written 
in terms of g-logarithmic function. A physical interpreta- 
tion has been given to q in the Richards' model. This has 
been possible due to the microscopic model presented in 
Ref. [21]. 
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